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Abstract. The CO-H2 conversion factor in galaxies is typically described as bi- 
modal: one value for discs and quiescent regions, and another (lower) value for mergers 
and starbursts. In this proceeding, I will describe both empirical observational evidence 
that the conversion factor varies with physical environment, as well as a theoretical 
model which aims to understand the physical processes which drive these variations. I 
present a functional form for Xco which can be applied to observations ranging in scale 
from ~ 70 pc to galaxy-wide scales, and show the consequences of the application of 
this model to the Kennicutt-Schmidt star formation law. 



1. Introduction 

Deriving an H2 gas mass from a giant molecular cloud (GMC) or galaxy full of clouds 
involves the usage of tracer molecules such as '^CO (hereafter, CO). This is because H2 
has no permanent dipole moment, and its first quadrupole moment lies ~ 500 K above 
ground, significantly warmer than the typical ISM temperature in a GMC. Converting 
from a CO emission line strength to an underlying H2 gas mass involves the usage of a 
CO-H2 conversion factor. 

The CO-H2 conversion factor goes by two names in the literature: XccQ and aco- 
The former has units of cm^^K-km s~^ while the latter M0 pc~^(K-km s~')~^ The 
two are equivalent, and simply related via Xco= 6.3 x lO^^aco in the aforementioned 
units. Since both Xco and aco are used in the literature, I will present the results of 
this proceeding in terms of both units. 

Observationally determining an X-factor requires deriving an independent mea- 
surement of H2 gas mass, and comparing this to an observed CO line flux. Via a variety 
of methods (some of these are summarised in the seminal conference proceeding by 
ISplomon & Barrett 1991), it has been found that the conversion factor is roughly con- 
stant in the Galaxy and nearby galaxies (when the metallicity is of order solar) with 
value Xco= 2 - 4 x lO^^cm-^/K-km s-\ or aco* 3 - 6 M© pc'^ (K-km s"') ' (for 



a re ference list of these observations, please see the introduction of Naray anan et al. 

SoTibh . 



However, despite the seeming constancy of the X-factor in relatively "normal" 
GMCs, observations of nearby ultialuminous infrared galaxies (ULIRGs) show that 
the usage of a Milky Way Xco in these environ ments would cause the infe rred H2 gas 
mass to exceed the measured dynamical mass ( Downes & Solomonlll998h . Hence, in 



'Xco is also referred to as the A'-f actor. I will use Xco or the X-factor arbitrarily in this proceeding. 
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ULIRGs, the X-factor must be lower than wit hin the Galaxy. In Figure [T] we show 
the distribution of Xqo and aco values from the'Pown es & SolomonI (Il998h survey, as 
well as a shaded region for the typical range of Galactic values. There is a dispersion 
amongst the ULIRG X-factors, though on average they are lower than the Galactic 
mean. Despite the observed dispersion in ULIRG conversion factors, the fact that they 
are lower on average than the Milky Way value has led the community to largely adopt 
a bimodal picture of Xco- a value of aco * 4 for disc galaxies, and aco ~ 0.8 for 
starbursts and mergers. 

The ramifications of a bimodal form of th e CO-H2 con version factor are signif- 
icant. As an example, iDaddi et ahl ( 2010a ) and iGenzel et all (HoiO) demonstrated the 
effects of assuming a bimodal Xqo on the Kennicutt-Schmidt (KS) star formation rate- 
gas density relation. While in Ssfr - ^^co space, galaxies at low and high-z lie on an 
arguably unimodal relation, the introduction of a bimodal Xqo to convert the CO lumi- 
nosity to an H2 gas mass causes the KS relation to become bimodal. That is, mergers 
and discs lie on different track on the KS relation when utilising a different Xqo for 
each. Th is is shown explicit ly in Figure El where I h ave compiled the galaxies pre- 
sented in iDaddi et all (1201 Oah and iGenzelet all (l2010h . and plotted (in the left panel) 
the ZsFR-2]co relation, and (in the middle panel), the Esfr-Sh2 relation which utilises 
a bimodal X-factor. The mergers (squares) and discs (circles) occupy different tracks. 
This has given rise to a terminology in the literature, as well as at conferences (including 
this one!), that mergers and discs have different "modes" of star formation. Mergers, 
according to this relation, form stars more efficiently (i.e. needing less gas to sustain 
a given SFR) than discs. We will neglect the right panel of Figure [2]for now, and will 
return to it shortly. 



2. The Case Against a Bimodal Conversion Factor: Empirical Observational Ev- 
idence 

The usage of a bimodal CO-H2 conversion factor raises a number of difficult questions. 
For example, if one ascribes a mean "merger/starburst" value of Xco to all mergers, 
should this value vary with the interaction stage of a merger? What about with close 
pairs which do not have overlapping discs? Does the same "merger/starburst" value 
apply for mergers of all gas fractions, orientation angles, star formation rates, mass 
ratios and galaxy masses? Does this merger value evolve with redshift, or should one 
utilise the locally-calibrated X-factor for mergers at high-z (which are presumably more 
gas rich and have higher star formation rates than their local counterparts)? Which 
value should one use for z ~ 2 discs, which are forming stars at rates comparable to 
local mergers? 

An alternative picture to the bimodal X-factor is one in which the CO-H2 con- 
version factor varies with the physical conditions in the ISM. Intuitively, this scenario 
makes sense. The fraction of hydrogen in molecular form, fraction of carbon in the 
form of CO, and escape of CO radiation from GMCs is all dependent on the physical 
conditions of the ISM. So, why shouldn't the CO-H2 conversion factor also depend on 
the physical conditions? In fact, there is substantial empirical empirical evidence that 
it does. 

First, as is shown in the Appendix of iTacconi et al.l (l2008h . X ro is seen to vary 



smoothly with the H2 surface density of galaxies. This was found by lOstriker & Shetty 
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Figur e 1 . His togram of the CO-H2 conversion factors measured for local ULIRGs 
by Do wnes & S olomon ( 1998) (solid black line). The yellow shaded region denotes 
the typical values found in the Milky Way. Despite a dispersion of conversion factor 
values for mergers, a singular value of aco=0.8 is typically assumed for merging 
galaxies. 



(|201lh to be well fit by a relation aco~ ^m^^ ~ ^co^^ where Zh2 is the surface 
density of molecular gas, and Wco is the observed CO intensity (in K-km s~^). 

Second, the Z-fac t or appears to vary with the metallici t y of the gas. Recent work 
bv lBolatto et~al] (l2008h . lLerov et allfcoil) and lGenzel eTall (I2OII1) has shown tha t Xm 
incre a ses with me tallicity with power X oc (0/H)~* where b = I - 2.7 ( Arimoto et al. 



Il996l : Ilsraeilll997h . The physical motivation behind this sort of scaling is that in low 
metallicity clouds, H2 can self-shield for survival from photodissociating radiation 
whereas CO requires some amount of dust to protect it. Hence, the fraction of CO- 
dark H2 mass rises as metallicity d ecreases IWolfire et al.ll201CH) . 

Right away, the results from ' Ostriker"& Shettvl (I201lh alongside the aforemen- 
tioned results regarding the variation of Xqo with metalhcity imply an empirical rela- 
tion which looks like: 



^co aco oc W^J^Z 



(1) 



3. A Theory for the CO-H2 Conversion Factor 
3.1. Numerical Models 

We will utilise numerical simulations to try to understand the origins of variations in 
the CO-H2 conversion f actor. The principle equations for this model are described in 
Narayanan et al. (1201 lal ). and I encourage the reader to refer to this paper. The descrip- 
tion here will be necessarily abbreviated owing to space limitations. We first model the 
evolution of idealised galaxies hydrodynamically utilising the publicly available SPH 
code GADGET-3 (Springcl 2005). We model isolated disc galaxies, and galaxy mergers 
over a range of masses, merger mass ratios, and redshifts. We smooth the physical 
conditions onto an adaptive mesh for the purposes of radiative transfer. 
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Figure 2. Kennicutt-Schmidt star formation relation in observed galaxies. Cir- 
cles and triangles are local discs or high-z BzK galaxies (inferred high-z discs), and 
squares are inferred mergers (local ULIRGs or hig h-z SMGs). Colours deno te sur- 
veys galaxies are pulled from and are described in iNaravanan et all (1201 Ibh . Left: 
SFR surface density vs. velocity-integrated CO intensity. In terms of the observ- 
able, Wco, a unimodal SFR relation exists. Centre: When applying an effectively 
bimodal Xco (cco = 4.5 for local discs, 3.6 for high-z discs, and 0.8 for mergers), 
the resulting SFR relation is bimodal, with mergers occupying a track of higher star 
formation efficiency (defined as SFR/Mh2)- The solid and dotted lines overplotted 
are the best fit tracks for each "mode" of star formation as published in lDaddi et al.l 
(1201 Oal) . Right: Resulting SFR relation when applying Equation|2]to the data. When 
Xco varies smoothly with galactic physical properties, the result is a unimodal SFR 
relation. In this case the power-law index in the relation is approximately 2 (when 
assuming that each galaxy has metallicity Z' = 1). 
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^The H2 fraction is calculated following the prescriptions of iKrumholz et all (12008 , 

20091) which balances the photodissocation by Lyman-Werner band photons against 



growth on dust grains. We assume that all H2 gas is assumed to be in spherical GMC. 
Because some cells (which contain little H2 mass) can be rather large in the adaptive 
mesh, we establish a floor surface density of 100 M0pc"^. In our model mergers, 
the bulk of the molecular mass is in GMCs which are above this floor value. We 
move beyond the 10,000 K floor temperature established in gadget-3 by calculating 
the temperature of the molecular gas as a balance between line cooling, heat ing by the 



grain photoelectric effect, cosmic rays, and energy exchange with dust (.Krumholz et al. 
20TTh . 



To calculate the CO emission line fluxes, we first determine the emergent ra- 
diation from each GMC utilis ing the publicly available escape probability code of 
iKrumholz & Thompson ( 2007h. We then ru n the 3D Monte Carlo radiative transfer 
package turtlebeach ( Narayanan et al.ll2b08h to calculate the radiative transfer through 
the galaxy. With modeled CO line fluxes, and knowledge of the H2 content in our model 
galaxies, we are prepared to study the CO-H2 conversion factor in our simulations. 



3.2. Results from Simulations 

At a given metallicity, Xqo decreases with increasing gas temperature and velocity 
dispersion. The reason for this is straightforward. The velocity integrated CO line in- 
tensity enters Xqo in the denominator. As gas temperature increases, the peak intensity 
of the CO emission line, which is in the Rayleigh- Jeans limit, also increases. Similarly, 
because GMCs are typically optically thick, as the gas velocity increases, the velocity 
integrated CO line intensity increases. 

This is the reason Xco is lower in galaxy mergers than the Galactic mean value 



(iNarayanan et al. 2 011a) . In a galaxy merger, the velocity dispersion of the model 
GMCs in our simulations increase from the virial velocity of the GMCs due to the tur- 
bulent nature of the gas in a galaxy merger. The velocity dispersion within a GMC 
in the merger can be as large as 50 - 100 km s~' (compared to a few km s~^ in a qui- 
escent cloud). Similarly, the gas temperature goes up in GMCs. While the typical 
temperature of a quiescent GMC is of order ~ 10 K, representative of the floor tem- 
perature established by a Galactic cosmic ray heating rate, the temperatures in GMCs 
in galaxy mergers can be of order - 50 - 100 K. The reason for this is due to en- 
ergy exchange with dust. When gas densities rise above ~ 10^ cm"^, dust and gas 
exchange energy extremel y efficiently, and the gas temperature rises to the dust tem- 
perature ( Goldsmiffil |200 ih . In a merger, not only are these large gas densities typical 



in GMCs, t he dust temperature is e levated due to the increased star formation rate (see 



Figure 2 of INarayanan et al.ll2011al) . A consequence of the increased gas temperatures 



and veloc ity dispersions in galaxy m ergers is a decreased Xqq. This may explain the 



results of Downes & SolomonI ([1998) and Figure [T] 

A second physical effect controlling the conversion factor is the gas phase metallic- 
ity. As explained before, in low-metallicity gas, the carbon resides predominantly in the 
form of CI or CII (rather than CO) due to photodissocation, whereas H2 can self-shield 
to survive. In this regime, the velocity-integrated CO intens ity decreases for a given Ht 
gas mass, and the X-facto r rises. MHD models of GMCs bv lGlover & Mac L ow (201 ij) 
and IShettv et all (l2011bUal) . as well as cosmological galaxy evolution simulations (R. 



Feldmann et al., submitted) find similar trends. 
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While our model isolated discs have mean Z-factors comparable to the Milky 
Way, and our model z=0 major mergers have X-factors consistent with the observed 
range presented in Figure [T] there is no "disc" value and "merger" value of Xqo- In 
a minor merger, for example, the rise in gas temperatures and velocity dispersions are 
not as extreme as they are in 1:1 major merger. Hence, the rise in velocity-integrated 
CO intensity is not as large, and the X-factor does not decrease as much. Similarly, 
in high-z disc galaxies, large gas clumps that are forming stars rapidly due to grav- 
itational instabilities in a gas-rich environment have larger temperatures and velocity 
dispersions than quiescent z=0 discs. Hence, they have X-factors lower than the locally 
calibrated "quiescent/dis c" value. There is no "merger" value and "disc" value of Xco 
(iNarayanan et alJIlOUbh . 



The good news is that the CO-H2 conversion factor can be observationally cali- 
brated. While the gas temperature and velocity dispersion are not easily observable, in 
our models, the CO surface brightness (line intensity) serves as a good proxy for the 
product of these quantities. That is, larger CO line intensities tend to correlate with 
larger velocity dispersions and temperatures. Fitting Xqo in our simulations against the 
CO line intensity, Wqo, and metallicity, we get: 

6.75 X 1020 X (Wco)-"-^' 
^co (2) 

10.7 X (Wco)"°'' 
"CO = (3) 

where (Wco) is the luminosity-weighted CO intensity over all GMCs in a galaxy, and 
Z' is the metallicity in units of solar. In the limit of a uniform distribution of luminosity 
from the ISM in a galaxy, (H^co) reduces to L^^/A, or the CO intensity (A is the area 
observed). Note the similarity of our model fit for Xqo, and what was derived from 
empirical observations in § |2l Hence, with a metallicity estimate and a CO surface 
brightness, one can calculate the X-factor observationally. I note that there is an implicit 
ceiling value of Xco - 4x lO^" cm^^/^-km s^^aco- 6.3 M© pc"^ (K-km s"')"') when 
Z' = 1 in Equations 2 and 3. 



4. The Ramifications of This Model Form for Xco 
4.1. The Kennicutt-Schmidt Star Formation Relation 

The usage of Equations 2 and 3 has an immediate consequence for the observed Kennicutt- 
Schmidt star formation relation. Recall that when a bimodal conversion factor is used 
for mergers and discs, a bimodal KS relation results with diff'erent normalisations for 
mergers and discs (centre panel of Figure El When utilising our continuous form for 
Xco, this bimodality goes away. 

We now turn to the right panel of Figure |2l where we hav e applied our mode l 



form for Xco to the observed data in the iDaddi et al.! ( 2bl0 al) andl Genzel et al.l ( 2010l) 



compilation. The bimodal relationship is reduced to a unimodal one with best fit index 
is a relationship Ssfr ~ ^mor ^^^^ picture, there are no different "modes" of star 
formation in mergers and discs. On average, mergers do form stars more efficiently 
than disc galaxies as they lie preferentially toward the high Xsfr regime of Figure |2] 
However, if one examines the central region of the fit, where both high-z discs (filled 
blue circles) and local galaxy mergers (open black squares) reside, it is evident that 
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Figure 3. Histograms of the derived conversion factors for all the galaxies in Figure|2] 



they have similar star formation efficiencies (defined as SFR/Mh2)- A key point of this 
model is that if the physical conditions in a disc and a merger are the same, they will 
have the same conversion factor 



4.2. The Typical Conversion Factors for Different Galaxy Populations 

This point is made more clearly in Figure [3l where we show the conversion factors 
we derive for the galaxies in Figure |2] utilising Equations [2l Here, we see that on 
average mergers have lower conversion factors than discs. The high-z SMGs (blue 
dashed line) have the lowest Xqo values, followed by low-z mergers (red dotted line). 
However, there is substantial overlap in the conversion factors, especially between the 
low-z mergers and high-z discs (black solid line). 

This model clearly predict, then, that attempts to determine conversion factors for 
high-z discs will evidence a lower average value than the G alactic mean value . A few at- 
tempts at constraining Xco in high-z BzK galaxies exist. iDaddi et al.l ( 2010bh developed 



dynamical models for two high-z BzK galaxies. After subtracting off the measured stel- 
lar and assumed dark matter masses, they were able to derive an Xco factor by relating 
the remaining (presumably) H2 mass to observed CO luminosity. This method recov- 
ered a mean X-facto r of ~ 2 x 10^" cm"^/K-km s~' , or aco ~ 3 pc~^ (K-km s"^)"^ 
Magdis et all (1201 ll) recovered similar results via dust to gas ratio arguments. Our 



model predicts that further measurements of the conversion factor from high-z discs 
will find a mean value a factor of a few lower than this (Figure O. 
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5. Conclusions 

The principal result from this work is that the CO-H2 conversion factor in galaxies 
varies with the physical conditions in the galaxy. It is not bimodal with the global 
morphology of the galaxy. Rather, the physical environment sets the CO-H2 conversion 
factor. 

The dominant drivers in the conversion factor are the gas temperature, velocity 
dispersion and metallicity. While the former two are difficult to observe, they can be 
parameterised by the observable CO intensity. With this, the CO-H2 conversion fac- 
tor can be described as function of CO intensity and gas phase metallicity. We show 
that this function can be derived both from empirical observational results, as well as 
numerical modeling. 
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